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Abstract. We study the possibility to measure the elastic = J ftp, '(/'(2S), '(/'(STTO), X2c) scattering cross sec- 

tion in the reaction p+d^'P+Usp and the elastic D{D)N scattering cross section in the reaction p+d D~ D'^psp. 
Our studies indicate that the elastic scattering cross sections can be determined for <f momenta about 4-6 GeV/c and 
_D/_D momenta 2 -5 GeV/c by selecting events with pt > 0.4 GeV/c for <f's and pt(psp) > 0.5 GeV/c for Z)/I)-meson 
production. 



PACS. 25.43. -l-t Antiproton-induced reactions - 
Decays of charmed mesons 



14.40.Lb Charmed mesons - 14.65. Dw Charmed quarks - 13.25.Ft 



1 Introduction 

Apart from the light flavor (qq) quark physics and their hadronic 
bound states the interest in hadronic states with strange flavors 
(ss) as well as their mutual interactions and properties in the 
medium has been rising continuously in line with the devel- 
opment of new experimental facilities [|l]]. In addition to the 
strange quark sector also the charm quark degrees of freedom 
have gained vivid interest especially in the context of a phase 
transition to the quark-gluon plasma (QGP) where charmed 
meson states should no longer be formed due to color screen- 
ing [§J3)0|. However, the suppression of J/tf' mesons in the 
high density phase of nucleus-nucleus collisions might also be 
attributed to inelastic comover scattering (cf. [H^^ and Refs. 
therein) provided that the corresponding J/tf'-hadron cross sec- 
tions are in the order of a few mb [||. Present theoretical esti- 
mates here differ by more than an order of magnitude espe- 
cially with respect to J/'f'-meson scattering. Also the J/^FN 
cross section is not known sufficiently well since the photopro- 
duction data suggest a value of 3^ mb [|o|, while the charmo- 
nium absorption on nucleons (at high relative momentum) in 
p + A and A + A reactions is conventionally fitted by 6-7 mb 
[0|n]] . In short: the present status of our knowledge on the in- 
teraction of charmed mesons with nucleons - especially at low 
relative momenta - is very unsatisfying. 

In this work we explore the perspectives of charmed meson 
- nucleon scattering in the pd reaction because i) antiproton 
annihilation will produce charmed mesons with rather low mo- 
menta in the laboratory and ii) the momentum distribution of 
the spectator nucleon in the deuteron is very well known. We 



mention that related experimental studies might be carried out 
at the future 'glue/charm factory' at GSI that is presently un- 
der discussion [p^]. An alternative possibility to measure the 
J /^N elastic scatternig cross section using pion beams in the 
reaction Tr+d ^ J/ipPP has been discussed recently by Brod- 
skyand Miller 

Our work is organized as follows: In Section 2 we will 
explore the perspectives for resonance production channels of 
J/tf', 1^/(25), 1^/(3770), X2c and their rescattering on the spec- 
tator nucleon including their decay to dileptons. In Section 3 
nonresonance production channels will be studied in the pd re- 
action with an emphasis on DD rescattering and the option to 
gate on charmonium momenta by triggering on an additional 
energetic photon or pion. Section 4 concludes this study with a 
summary. 



2 Resonance production and rescattering in 
the reaction p+d^<P+n 

We here examine the possibility to measure the elastic iPN 
(<1> = J/?/;, ?/'(25), ■(/'(3770), X2c) scattering cross section in 
the reaction 

p+d^<P+nsp^X+nsp, (1) 

where <1> is produced by resonance fusion pp ^ (p and X the 
decay product of the resonance. The main contribution to the 
total cross section of the reaction ([l|) comes from the specta- 
tor term (see diagram a) in Fig 1.), which is dominant if the 
momentum of the proton-spectator from the deuteron is below 
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100-150 MeV/c. The amplitude corresponding to the diagram 
a) of Fig. 1 can be written in the deuteron rest frame as 



'Mr), (2) 



where f{pp <P X) is the Breit-Wigner amplitude, p2 the 
proton momentum and !fd(r) the deuteron wave function. 



Table 1. Charachteristics of the transitions pp ^ $m vacuum and for 
a deuteron and carbon target. 







^{2S) 


*'(3770) 


X2c(lp) 


M^{MeV) 


3097 


3686 


3770 


3556 


r4.{MeV) 


0.087 


0.277 


23.6 


2 


kcn.{pp)iMeV/c) 


1232 


1586 


1634 


1510 


Tiab(MeV) 


3250 


5400 


5700 


4845 


Br{pp) X 10"" 


21.4 


1.9 


2 - 0.2 (?) 


1.0 


a{pp^$)(pb) 


5.15 


0.28 


0.7-0.07 


0.27 


Br{e+e-) x 10"* 


600 


85 


0.112 




Br{'yJ/tp) 








0.135 


S.F.{d) X 10"^ 


1.07 


2.6 


221 


19.2 


S.F.{C) X 10"^ 


0.56 


1.24 


113 


11.2 



All resonances J/ip, ip{2S),ip{377Q) and X2c are quite nar- 
row and the decay length for each of them produced in the re- 
action (|l]) is much larger than the average distance between 
two nucleons in a deuteron (see Table 1). Thus the produced 
meson can rescatter elastically on the neutron spectator, trans- 
fering momenta larger than a few hundred MeV/c. In this case 
the spectator term will become very small and the contribution 
from the rescattering term become dominant. 

Each resonance (produced in the pd reaction) has a rather 
high momentum with respect to the spectator nucleon such that 
we can describe the rescattering amplitude (see diagram b) in 
Fig. 1) within the framework of the eikonal approximation (cf. 
Ref. [|lll||l6|]) as 

Mb = -fipp ^^^X)Jdh e-'P^''0{z)r{h)Wd{r), 

(3) 

where the z- axis is directed along the ^-momentum. The elas- 
tic 'Pn - scattering amplitude is related to the profile function 
r{h) by the standard expression 



/(<??! ^ ^n) = / d^h e-'^i^rih) 

znik J 



(4) 



with k denoting the laboratory momentum of the <?-meson and 
q the momentum transfer. 

The probability to detect the <P through the decay channel 
X at solid angle df2d in coincidence with the spectator of mo- 
mentum p2 is related to the differential cross section as 



d^a{pd ^ <Pn ^ Xn) 
df2d^P2 



= \Ma + Mb 



(5) 



In order to simulate events for the reaction (|l|) we use the 
Multiple Scattering Monte Carlo (MSMC) approach. An earUer 



X 



a) 



b) 



Fig. 1. The diagramms for ^-meson production in the reaction (|l|) 
without (a) and with ^ A'^ rescattering (b). 



version of this approach - denoted as Intra-Nuclear Cascade 
(INC) model - has been applied to the analysis of 77 and oj 
production in pA and pA interactions in Refs. [ p^|Jl8| ] . Recently 
this version of the INC model has been extended to incorporate 
in-medium modifications of the mesons produced [ [l9po| ] in 
hadron-nuclear collisions. 

However, the INC model is valid only for medium and heavy 
nuclei and cannot be directly applied to deuterons. In order to 
perform simulations of scattering events in the case of antipro- 
ton - deuteron interactions we use the Monte-Carlo approach 
for the single and double scattering terms, which are propo- 
tional to \Ma\'^ and |Aff,p of Eq. (q]), respectively. Note that 
the interference between the spectator and rescattering ampli- 
tudes is only important in a narrow region of spectator mo- 
menta, where both contributions are of the same order of mag- 
nitude. For a more detailed discussion of this point see, e.g., 
Refs. [[l|,|l5|jl|,p. 

The probability for the produced meson to rescatter on the 
spectator nucleon then can be found in a standard way by 



with 



(6) 



(7) 



Necessary parameters for a MC simulation of rescattering are 
the elastic scattering cross sections and slope parameters 
h for the differential cross sections da/dt, which are approxi- 
mated by 

da/dt = Aexp{bt), (8) 

where t is the momentum transfer squared. These parameters 
as well as the masses of the rescattered particles determine the 
momentum and angular distributions of the particles in the final 
state. 

An important point for resonance production in pd reac- 
tions is the rather strong dependence of the Breit-Wigner cross 
section \f{pp (P X)\'^ on the initial proton momentum 
which leads to a rather strong suppression of the cross section 
on the deuteron or on heavier targets. The suppression factor in 
the maximum of the cross section is of the order (see Ref. [^) 

S.F. ~ nr^mpf /{kpAI^), 

where F^,, M4, denote the vacuum width and mass of the pro- 
duced meson, whereas fc^? is the target Fermi momentum. More- 
over, after Fermi smearing over the target momentum distribu- 
tion the Breit-Wigner production cross section becomes much 
wider and asymmetric (see e.g. Ref. [^^). 
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All the effects of Fermi smearing have been taken into ac- 
count in the MSMC calculations explicitly by convoluting the 
resonance production cross section with the proton momentum 
distribution in the deuteron (or heavier nuclei). In a deuteron 
this distribution is given by the Fourier transformed deuteron 
wave function squared while in medium and heavy 

nuclei the local density approximation for the Fermi distribu- 
tion is used. 
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Fig. 2. Fermi smearing of (j(pp ^) for (p = J/tp,tp{2S) and 
'!/j(3770) in deuteron (middle part) and carbon targets (lower part). 
The upper part shows the elementary cross sections on a hydrogen 
target. 
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Fig. 3. Distributions of J ftp mesons from the reaction ( |lC| ) in longi- 
tudinal (upper left) and transverse (upper right) momentum for Tiab = 
3.25 GeV. The distribution of the decay products e^/e^ in the total 
momentum is shown additionally (lower left) as well as in the po- 
lar angle 6 (lower right). The solid and dashed curves describe the 
contributions of the spectator and rescattering terms, respectively. The 
J/ipn elastic scattering cross section was assumed to be 1 mb and the 
differential slope parameter was taken as b ~1 (GeV/c)~^. 



In Fig. H we present the cross sections a{pp J/tp)^ 
(j{pp ■ip{2S)) and a{pp — > -0(3770)) as a function of the 
antiproton kinetic energy for a hydrogen target (upper part), 
for a deuteron target (middle part) and carbon target (lower 
part). In Fig. 2 two excitation curves for J/ip production on 
the deuteron are shown to demonstrate the dependence on the 
choice of the deuteron wave function: the dashed line corre- 
sponds to the Paris model while the solid curve is calcu- 
lated using the Hiilthen wave function 

Mr) = N (e-"'- - e-^^) /r, (9) 

where N is the normalization factor, a = ^me ~ 46.5 MeV/ c, 
P ~ 5.2a. One can see that this dependence is quite small 
since the two results practically overlap. 
Let us now consider the reaction 

pd — *■ J/ipn e~n (10) 

In Fig. ^ we show the distributions of J/ip in longitudinal (up- 
per left) and transverse (upper right) momentum as well as the 
momentum distribution of the decay products e+/e^ (lower 
left) and their distribution in the polar angle 6 (lower right) 
for Tiab = 3.25 GeV. The solid and dashed curves describe 
the contributions of the spectator and rescattering terms, re- 
spectively. The J/ij] n elastic scattering cross section here was 
assumed to be 1 mb while the slope parameter was taken as 
6 ~1 (GeV/cV^, which is almost the same as in the reaction 
7P — > J/ tpp [gsj] at roughly the same s ~ 20 GeV^. 



It is seen that after rescattering on the neutron the J/^/; 
looses longitudinal momentum (upper left part), but gains trans- 
verse momentum (upper right). The momentum spectrum of 
the leptons from J/t/j decay is completely flat (lower left part), 
while their angular distribution (lower right) has a broad maxi- 
mum around 25°. 

It is clear that performing a cut for pt{J / ii) > 400 MeV/c 
or pl{J/4>) < 3.6 GeV/c one can determine essentially the 
overall magnitude of the J/tpn elastic cross section. This is 
demonstrated in Fig. 4 where we show the number of events as 
a function of the cut in the J/ip transverse momentum. The 
dotted and dashed curves describe the tails of the spectator 
momentum distribution for the Paris and Hiilthen models of 
the deuteron wave function, respectively. For pt{J /ili) > 400 
MeV/c both curves drop far below the solid histograms which 
describe the contribution of the J/tpn rescattering term. The 
relative number of events with rescattering is propotional to 
(TeiiJ/ipn) and for the cut pf( J/ t/;) > 400 Me V/c is about 0.1 
% using (Jei[J /ipn) = 1 mb. 

In a similar way we can evaluate the resonance production 
and rescattering of particles from the ^ family as ■!/'(2S'), -0(3770) 
and X2c(lp) which can be detected via their decays l^l^ , DD 
or "fJ/tl), respectively. The suppression factors due to Fermi 
smearing of the cross sections ^ in deuteron and carbon 
targets for all particles are shown in Table 1 . 

The case of the ip{2S) looks much less promising than the 
case for J/ip, because its signal in the l^l^ channel will be 
smaller by almost two orders of magnitude. A much stronger 
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Fig. 4. The number of events for reaction ( [lo| ) as a function of the pt 
cut for pt{J /ijj) > Pf"*. The dotted and dashed curves describe the 
contributions of the spectator term for the Paris and Hiilthen models 
of the deuteron wave function, respectively. The solid curve describes 
the rescattering term for Oci{J /^}jn) = 1 mb and a slope parameter for 
the differential cross section 6=1 (GeV/c)~^. 

signal of the ^{2S) will be in the channel J/i/itt+tt^, but in 
this case a measurement of the '0(25') transverse momentum is 
difficult. 

The signal of the X2c(lp) in the channel jJ/ip will be al- 
most the same as the signal of the J/^ in the dilepton chan- 
nel. Despite of its smaller production cross section in pp the 
X2c(lp) is not so strongly suppressed by Fermi smearing as in 
case of the J/ip (cf. Table 1). Therefore, if the photon from the 
X2c(lp) decay can be detected in coincidence with the dilep- 
tons from the J/tp decay, a measurement of the X2c{^p)n elas- 
tic cross section can also be performed by selecting events with 
Ptix) > 400MeV/c. 

It is, furthermore, interesting to consider the reaction 



pd ^p{3770)N ^ DDN. 



(11) 



Up to now the resonance ^'(3770) has been seen only in the 
charge neutral mode in e~^e~- collisions. Therefore, its isospin 
is not yet known and a detection of DD production in pn would 
unambiguously fix its isospin as 1 . If the DD pair will be pro- 
duced only on a proton, then its isospin will be 0. 

The coupling of the V^(3770) to the pp channel is not known. 
Having in mind that close to the ■0(3770) the resonances X2c(lp) 
and tp(2S) have a branching ratio BR{pp) ~ 2.10"^ (despite 
of their full widths being different by almost an order of magni- 
tude), it seems to be justified to assume i3i?('i/;(3770) pp) — 
2.10"''^ — 2.10~^. Then the cross section of its production on a 
nucleon in the deuteron will be about 15-150 nb. 

In Fig.BI we present the simulated momentum distributions 
of 0(3770) and D/D mesons for a{pp i/'(3770)) =150 nb, 
ct(V'(3770)7V -0(3770)^^) =10 mb and a slope parameter 
b{ip{3770)N)= 2 (GeV/c)-2. The solid and dashed curves de- 
scribe the contributions of the spectator and rescattering terms, 
respectively. The 07(3770) distributions in the longitudinal and 
transverse momentum are shown in the upper left and upper 
right part, respectively. The rescattering term dominates for 
PL < 6 GeV/c and pt >0.4 GeV/c. The pi and pt spectra of 
D/D mesons are presented in the lower left and lower right 
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Fig. 5. 



Momentum distributions in pl (left) and pt (right) for 
■0(3770) (upper) and D/D (lower) from the reaction (|l l|). The 
solid and dashed curves describe the contributions of the spectator 
and rescattering terms, respectively, with <j{pp tp{3770)) =150 
nb, a{tp{3770)N ^ V(3770)iV) =10 mb and a slope parameter 
b{Tp{3770)N= 2 (GeV/c)-2. 



part of the figure. The change of the longitudinal and transverse 
momentum spectra for ■0(3770) due to rescattering has also 
some influence on the pi and pt spectra of the D /D mesons 
especially for the lowest pi and the highest pt- However, the 
signal of rescattering is much cleaner in the pt distribution of 
the 0(3770) for pt >0.4 GeV/c: practically there is no back- 
ground from the spectator mechanism and the relative number 
of the events is 0.73 % in case of the cross sections quoted 
above. 



3 Nonresonance production and 
rescattering of J/ip and D/D 

Since there is a quite large suppression in the resonance pro- 
duction for narrow quarkonium states on nuclear targets, their 
nonresonance production (where a suppression factor is absent) 
might well be comparable to the resonant formation. In this 
Section we thus consider the rescattering of J/t/j and D /D in 
the deuteron where the mesons are produced in the nonresonant 
reactions 
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Fig. 6. Longitudinal and transverse momentum distributions of J/tp 
(lower part) and the proton spectator (upper part) from the reaction 
Am at Tiab =5 GeV. The solid curves are the contributions from the 
spectator term. The contributions from pion and J/tp rescattering on 
the proton spectator are shown by dotted and dashed curves, respec- 
tively. The parameters of the J/tpN scattering amplitude were taken 
the same as in Fig. 3. 



pd J/lp TT p 

pd ^ J 1'^ jn 
pd -> DDN. 



(12) 
(13) 
(14) 

The cross section of the reaction ( [l4l ) has been calculated 
in Ref. [|6| within the framework of the Quark-Gluon String 
Model. It grows from threshold up to piab — 7.5 GeV/c, where 
it reaches its maximal value of about 40 nb, and then decreases 
rather fast with energy as 

52[a(0)i;,-l]^ 

where a{0)s^ is the intercept of the Sc Regge trajectory 

a{t)s, ~ -1.82 + 0.5^. 

The cross sections of the reactions (|l2|) and ( p3| ) in the con- 
tinuum are not known. We assume that a{pd J/ip tt^p) is 
comparable with a{pd DDN) and take for our estimates 
30 nb at Tiab = 5 GeV. 
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Fig. 7. Correlation between the longitudinal momenta of J ftp's and 
photons at piab = 6, 9 and 12 GeV/c (l.h.s.) and the angular distribu- 
tions of photons at the same momenta in the reaction (O). 



The cross section for the reaction pp — > J ftp 7 can be esti- 
mated at Tiab = 4.845 GeV where it is determined through the 
resonance X2c{^p) (cf. Table 1); it is about 30-40 nb. Of course, 
in the continuum it is expected to be smaller. 

The longitudinal and transverse momentum spectra of J/ip 
and the proton spectator from the reaction (|l2|) at Tiab =5 GeV 
are shown in Fig. ^. The solid curves describe the contributions 
from the spectator term, the dottted and dashed histograms de- 
scribe the contributions from pion and J/ tp rescattering on the 
proton spectator The parameters of the J/tpN scattering am- 
plitude were taken the same as in Fig. 3. The initial antiproton 
energy was chosen to be above the J/tpn threshold, but slightly 
below the J/ tpinr threshold. Despite the kinematics of the re- 
action ( p^ is less resctrictive and the J/ ip has a broader distri- 
bution in transverse momentum, it is still possible - using the 
cut pt{J /ip > 0.7 GeV/c - to obtain a rather clean signal from 
the J/t/jN rescattering events. 

Similar considerations can be made for the reaction (pj|). In 
this case it is possible to use the photon momentum as an addi- 
tional trigger for J/t/m rescattering events, which can be used 
to fix the momentum of the J/t/j. The correlation between the 
longitudinal momenta of J/t/ys and photons is presented for 
different momenta of the antiproton in Fig. ^ (l.h.s.). As com- 
pared to the reaction ( [1^ ) one looses about 2 orders of mag- 
nitude in event rate due to the electromagnetic vertex, but the 
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PJGeV/c] P,[GeV/c] 




PJGeV/c] P,[GeV/c] 

Fig. 8. Longitudinal and transverse momentum distributions of D"- 
mesons (upper part) and proton spectators (lower part) from the reac- 
tion ( p^ at piab =7.5 GeV/c. The solid curves are the contributions 
from the spectator term while the contributions from D^^N rescatter- 
ing on the proton spectator are shown by the dashed histograms. The 
elastic D"p scattering cross section and the slope parameter b were 
assumed to be 10 mb and 2 (GeV/c)~^, respectively. 



descrimination from background events is rather clean. The an- 
gular distributions of photons in this reaction are also presented 
in Fig. 7 (r.h.s.) and show that the photons are predominantly 
emitted to forward angles. 

Let us now discuss the possibility to measure the DN and 
DN elastic scattering cross sections in the reaction (^4|). In 
Fig. 8 we show the longitudinal and transverse momentum dis- 
tributions of _D"- mesons (upper part) and proton spectators 
(lower part) at piab =7.5 GeV/c. The solid curves are the con- 
tributions from the spectator term while the contributions from 
D'^N rescattering on the proton spectator are shown by the 
dashed histograms. The elastic D'^p scattering cross section 
and the slope parameter b were assumed to be 10 mb and 2 
(GeV/c)~^, respectively. In the case of DN scattering we have 
used (Tei = 5 mb and the same slope parameter b. In this case the 
contribution of the D^p scattering term is smaller by a factor 
of 2. It is seen that the spectator term dominates for p^p < 400 
MeV/c. However, at psp >500 MeV/c the main contribution to 
the spectrum comes form the rescattering term. The individual 
contributions from D'^p and D~p rescattering can be separated 




40 80 120 160 200 240 280 320 360 

t [deg.] 

Fig. 9. The distributions in the azimuthal angle (p between the planes 
piab- Pd" and piab-Psp (upper part) and between the planes piab- 
Pq- and piab-Psp (lower part) for pt {pap) > 500 MeV/c in the reac- 
tion pd DDN at Piab = 7.5 GeV/c. 

using the correlation between the azimuthal angles of the two 
scattering planes put- Po/Qand piab-Psp- 

This correlation is presented in Fig. ^ for pt{psp) > 500 
MeV/c where we show the distribution of events in the az- 
imuthal angle between the planes piab- Pd" and piab-Psp in 
the upper part and between the planes piab- Pd- ™d Piab-Psp 
in the lower part. The shaded area describes the background 
from D^p (upper part) and D'^p rescattering (lower part). 



4 Summary 

In this study we have explored the perspectives of measuring 
the elastic cross section of charmed mesons with the spectator 
nucleon in resonant and nonresonant pd reactions. Our analy- 
sis within the MSMC approach indicates that the elastic scat- 
tering cross sections can be determined for <!>{= J/ip, tp{2S), 
7/1(3770), X2c) momenta about 4-6 GeV/c and D/D momenta 
of 2 - 5 GeV/c by selecting events with pt > 0.4 GeV/c for <?'s 
and Pt {psp) > 0.5 GeV/c for D / D-meson production. 

We mention that the inelastic cross sections of charmed 
mesons may be studied in pA reactions as analysed in Ref . [|7| . 
This opens interestin g pe rspectives for a future high energy an- 
tiproton storage ring [ |l2| ] . 
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